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The starting formula of Bosvieux and Friedeli for the force on an ion in a metal due to 
an applied voltage is shown to lead to the same description as the linear-response approach 
used in the field since its introduction by Kumar and Sorbello^. By this electromigration 
theory has become a unified theory. It follows after accounting for a treacherous trap term, 
which at first sight seems to be zero. Up to now Bosvieux and Friedel claimed to predict 
a completely screened direct force, which means that only a wind force would be operative. 
In addition the amount of screening has been calculated up to second order in the potential 
of the migrating impurity, using a finite temperature version of the screening term derived 
by Shami For a proton in a metal modeled as a jellium the screening appears to be about 
15%, which is neither negligible nor reconcilable with the old full-screening point of view. 

Introduction 

Electromigration is the motion of an ion in a metal under the influence of an applied voltage. Depending on the 
metal and the ion its motion is either to the anode or to the cathode. The total driving force on such an ion is known 
to be the result of two contributions, a direct force and a wind force 

F = Fdirect + F w i n d = (Zdirect + Z w i n d)eE = Z*eFj, (1) 

in which Z* is called the effective valence, which is a measurable quantity. The direct force is due to the direct action 
of the applied field E on the charge of the ion. The wind force comes from the scattering of the current carrying 
electrons off the ioni4 

Until 1962 there was a common belief that Zdirect was equal to the bare valence Zi of the ion. At most a small 
deviation from that value could arise from the electrons in the metallic environment of the ion. In 1962 Bosvieux 
and Friedel predicted a complete cancellation of the direct force due to screening effects, so that only the wind force 
would remain. By that prediction a controversy was born. It was not easy to decide matters by a measurement, and 
a satisfactory theoretical answer was lacking as well. 

In 1975 Kumar and Sorbello^ published an exact linear-response expression for the driving force. After that it was 
considered as being just a matter of a careful evaluation of that expression in order to settle the problem. It took 
quite some time to do so though. Finally Sorbello^, starting from the complicated treatment given by Rimbey et al&, 
predicted a screening of 10-30%, depending on the potential used for the ion. Support of this result has been given 
by the present author, starting from a much simpler description. 7 However, in the meantime new support was given 
for the full-screening point of view£ 

Interestingly, the original paper by Bosvieux and Friedel (BF) is referred to as being controversial on the one hand, 
while on the other hand their starting formula is referred to as the first quantum mechanical equation for the wind 
force, which is seen as a pioneering contribution to the field. Therefore we undertook an new evaluation of their 
starting formula, following the authors as closely as possible, up to a point where we came across a trap. By treating 
this trap properly the final evaluation ends up at the standard linear-response expression for the driving force. As far 
as the author knows Sorbello was the first who showed that such traps occur in the theory of electromigration. 9 

The screening itself has been calculated by Sorbellc-, who used one type of expression for it. Therefore we undertook 
the evaluation of Sham's second order expression written in a new form, the more so as a more recent result for the 
screening, ranging from to 100%, is rather inconclusive. 10 Sham's expression has never been evaluated numerically, 
and in his paper he just gives an order of magnitude comparison with his wind force expression. 3 This has led him to 
the conclusion that the screening is negligible. We find a screening of 5-25%, which is in agreement with Sorbello's 
results. 

Rydberg atomic units are used, in which the energy is in Rydberg, the distance is in Bohr (1 Bohr w 0.5 A), h = 1, 
and the electronic mass is i. 
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Basics of linear-response theory 

We first give the standard linear-response expression for the driving force on an ion at position R^r^ 

F = Z ie E - ieE v dte- at Tr{p(H) F op (t), }, (2) 

in which the cartesian label v runs from 1 to 3. The infinitisimally positive number a represents the adiabatical 
switch-on of the electric field, the operator p(H) is the grandcanonical density depending on the system Hamiltonian 
H, the force operator F op stands for 

F op ee -V Ri F = - VnA*3 ( 3 ) 

j 3 

and its time dependence refers to the Heisenberg representation 

F op (i) ee e iHt F op e- iHt . (4) 

The expression published by Kumar and Sorbello^ follows simply after a partial integration in Eq. (J2J) with respect 
to the time. 

All important studies of the driving force have been done for the electron-impurity system, the Hamiltonian H of 
which can be written as a sum of single particle Hamiltonians h, so 

H = J2h j with h = h +v = ho + ^V, (5) 

j a 

in which the summation in the electron-impurity potential runs over the positions R Q of the impurities. This allows 
for a reduction of the many body expression ijSJ) to the following single particle expression, 

¥ = Z i eE-ieE v J dt e~ at tr{n{h)[f\t), /]}, (6) 

where n{h) is the Fermi-Dirac distribution function in operator form 

The treatment of Bosvieux and Friedel (BF) 

We start to follow BF's treatment, by writing down their system Hamiltonian 

- A 

H=-J2^t + n^---^---^)- (8) 

In applying a small perturbation denoted as SV they write down the perturbed wave function ty. In their appendix 
BF employ the idea of switching on the field adiabatically. For the sake of clarity we give the explicit form of the 
perturbing potential due to an applied field, 

SV(t) =eEe Qt - (Y^rj -J2 z ^a) = SVe at = (5V (r) + 5V (R)) e at , (9) 

j a 

which is zero in the limit t —> — oo. The position dependent potentials have been defined for later use. BF's force 
expression reads as 

4> = - < V\Vi(V + 8V)\V >, (10) 

in which |^ > represents the state of the system as it develops from its unperturbed ground state 1-00 > due to the 
perturbation SV. In order to find \ty > we solve the time dependent Schrodinger equation 

i^=W(t)*(t) (11) 
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for the total Hamiltonian 

H(t) = H + 6V{t) (12) 
by working in the interaction representation for ^(t), defined as 

ttj(t) = e lHt ^{t). (13) 

The equation for ^i(t) becomes 

= e iHt(_ R + i^)m(t) = e mt SV{t)e- im ^i{t). (14) 
After integrating this equation one finds for ^(t) linearly in 5V 

m{t) = -ie- lHt f *'e^ t '(5t/(t')e^ ,;fft '*/(-oo) +e- iHt * 7 (-c3o). (15) 



With 6V(t) = 5Ve at , applying the substitution t — t' = s, and considering an arbitrary time in the present, so t = 0, 
this becomes 

/>oo 

tf(0) = * = -*/ dte- {lH+a)t 8Ve lHt ^> I {-<x>) + ^ I {-<x>). (16) 



If one calculates matrix elements with this 1^ > the factor e - lE o°° [ n the state |4 r /(— oo) > drops out so that just the 
ground state \ipo > remains. BF's Eq. (1.2) for |^ > is reproduced by inserting the complete set of eigenstates of the 
system Hamiltonian H, denoted by \tp n >. The only difference is the presence of the infinitesimal number a. BF have 
omitted it and just denote that in the summation over n the value n — 0, corresponding to the ground state, has to 
be excluded. It will become clear that this difference has a dramatic influence on the final results. 
Now we evaluate BF's force expression, Eq. (ll(Jf) . linearly in the applied field SV. 

/>oo 

$=-< o |V 4 ((5^)|Vo > +» / dt e~ at < ipo\(ViV)e- iHt 6Ve im \i/)o > + cc. 

Jo 

poo 

= Z,eE + i dt e~ at < Vol (V^), e - im SVe iHt >, (17) 
Jo L J 

in which we used Eq. Q16[l for > and Eq. for evaluating the matrix element < ipo\V i(5V)\ipo >. 

In the further evaluation of this equation we want to follow BF's evaluation given in their §1. The integral over the 
time in the first line of Eq. 117|l is carried out after inserting a complete set \ip n > of eigenstates. Further, use is 
made of the equality (VjV) = [V<, H]. This way one obtains 

$ = Z,eE J2(E n E ) < W„ > < J nl6 ^ l ° > + cc. 

ho — E n + ia 

= Z 4 eE - V < >< V„|<Wo > En ~ E °~ ia + ia + cc. 

^ ' hrx fr,__ 4- in 



Eq — E n + ia 
^— ' Eq — E n +ia 



= -ia ) < V>o Vi Vn > p ; ■ + cc. (18) 

^ Eq — E n + ia 

We want to comment on Eq. 1|18|) in view of the results of BF. First, in evaluating the matrix element 
< f/>o|[Vi, H]\ip n > in the first line of Eq. 1|18|) BF create in addition surface integral terms corresponding to Green's 
theorem, by that not appreciating the hermitian property of H . These surface integral terms represent flow of prob- 
ability out of the system, which is obviously zero for a finite system and for an isolated metal, as BF admit. But this 
flow is zero also for a metal carrying a steady electric current, which is denied in practice by BF. Secondly, in their 
treatment the last line of Eq. i|18[) is missing, because BF have put a = from the beginning. We want to point out 
that regarding this last line a treacherous, hidden trap in the formalism is involved, which has shown up earlier in 
electromigration theory^.. This term seems to be zero because of the proportionality with the infinitisimal number a, 
but we will show that this term in fact is a rich one. 
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In rewriting the last line of Eq. (|18f) we first write the energy denominator as coming from a time integral as follows 

/•OC 

$ =-a dt e- at Y] < M^i^n >< ipn\e- im 8Ve iHt \ip > + c.c. 
Jo 

= -a f dt e~ at <^ |[v 4 ,e- 4Ht ((5t/(r)+(5T/(R))e lfft ]|Vo > (19) 



Because H commutes with 5V(R) the second term reduces to ZieE, and one finds that 

/*°° 

$= Z ie E -a dt e~ at < ip \ e tHt \7 ie - lHt , 5V(r) \ipo > ■ 



(20) 



After a partial integration with respect to time, using that [Vj, H] = (ViV), and substituting SV(r) according to Eq. 
©, one finds 



= Zn E + / / dt e~ at < Vol e tHt (y r V)e- lH \ eE • ^ r A |^ > 



(21) 



We are back at the second line of Eq. I|17[l . but only the electron coordinates in 5V have survived the operations. 
The influence of the electric current on the total force is taken care of by the second term in the right hand side of 
Eq. I|21l) . In retrospection this result follows rightaway from Eq. H17|) as well, but we wanted to go along with BF's 
way of evaluation first. 

Eq. I|21|l is a very interesting result. It is precisely the zero temperature equivalent of Eq. @. This becomes even 
more clear if we write down the form which shows up after the reduction of Eq. (|21|l to single particle states denoted 
by \q >■ 



ZieE 



dt 



~ at zZ 

q 



<Q\ 



ill I 



ill t 



eE- 



\q>. 



(22) 



At T = the sum over the single particle states has a sharp cut-off at e q — ep. The finite temperature equivalent of 
Eq. (|2*2l can be written as 



Z t eE + i dt e~ at n ( € q) < l\ \e M {\7 lV )e- M , eE ■ r 



Z,eE -i dt e- at tr\ n{h) P(t),eE ■ r } = (Z, + Z wmd + Z scr )eE, 



(23) 



in which the Fermi-Dirac distribution n(e) has been inserted, see Eq. (JJJ). The force operator f* is defined in Eq. 
©. Clearly, Eq. 123|l is completely equivalent to Eq. © of the present text. By this electromigration theory can be 
considered as being unified. Apparently, BF's starting formula was correct, but these authors did not recognize its 
precise contents. This became even more clear recently^ BF used Eq. I|10|l for their result regarding the direct force 
only. They treated the wind force separately, applying a semi-classical standard approach in describing the current 
carrying electrons and accounting for the scattering of the electrons by the migrating impurity quantummechanically. 13 
The formal Eqs. J2J or l|23|) . or forms which have been shown to be equivalent to them, have been used since 
Kumar and Sorbello published their linear-response approach to the electromigration problem 2 . The result of all this 
research is, that the second term of these equations contains two contributions, as it has been indicated explicitly in 
the right hand side of Eq. (|23|) . One can be identified as the wind force. The other one implies some screening of 
the bare direct force ZieE. The wind force expression has been studied thoroughly and has been applied to calculate 
Zwind for both interstitial and substitutional impurities in many metals^^**^ Because the present paper is devoted 
to the settlement of the controversy regarding the direct force, we want to add a calculation of Z scr starting from 
Sham's contribution. 3 But first one more comment on Eq. Q23[l and its interpretation. BF claim to have proven that 
for an isolated system the force on an ion is zero^ The present author agrees that this force is zero, but this fact 
does not follow from an explicit proof, but from the knowledge that internally any applied field E is screened out by 
an electronic surface charge that is built up. This is a general result from the theory of electromagnetism. That is 
why then </>, which is proportional to E, is zero, because E = 0. 
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Calculation of Z acr 

If one evaluates Eq. (IL'MI) to lowest order in the potential v of the ion, for a jellium model of the metal, one finds 
thafeli 



= -— V(fc 2 -k-k') 

3m ^— ' 



9m 



n fc - rife/ 



(24) 



in which a fc label refers to a plane wave. It has been shown that the T — > limit is equivalent with Sham's 
expression. 17 Sham gave an order of magnitude estimate by comparing it with the wind force expression. Z win( i is 
proportional to the transport relaxation time r of the system, while Z scr is proportional to the inverse of an energy, 
for which the Fermi energy can be chosen. By that he came out at a ratio of (cft) -1 sa 0.01, being negligible. 
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FIG. 1: The amount of screening represented by Z scr according to Eq. [I2.H . with a = r , for the screened Coulomb potential 
and for two square well potentials. 

A numerical evaluation of Z scr becomes possible if one employs the spherical wave expansion for a plane wave, 
converts the summations over the wave vectors to integrals and carries out the angular integrals over the directions 
of the wave vectors. After using the relation between k 2 and the energy e/- one ends up at 



3TT 2 m J 



dei 



dn k _ n k -n k , 



in which the function fe(k, k') is defined as 

fi(k, k') = e ky /e^vi(k' k) 



(21 + l)kv e (k', k) - 2(1 + l)k'v e+1 (k', k) 



containing the information about the ion potential through 



vt(k', k) 



r 2 dr ji(k' r)v(r)jf (kr) 



(25) 



(26) 



(27) 



The integrand has to be treated with care when = because then the denominator attains the value a 2 which 
would imply 'singular' behaviour. However, precisely then the numerator becomes zero, because ]im e ,—> ek (nk — 
«fc')/( £ fc — e k') — * §t^- The crucial part of the integrand lies in the square around the point (eu^ey) — (eF,ep). In 
studying the Z scr integral it appears that in that square one has to keep the Fermi-Dirac distribution function in its 
finite temperature form. A calculation for o « r 1 has been presented elsewhere^ Here we want to follow Sham 
more closely. He replaced a by r _1 , by that accounting for the entire distribution of impurities in the metal and for 
the possible presence of other scattering mechanisms. 



6 



The result of a numerical evaluation for different ion potentials is shown in FIG. ^ I n addition to a screened 
Coulomb potential^ 

nr 1 -Xr 

v(r) = (28) 

with the Thomas-Fermi screening parameter A, square well potentials were employed in the same spirit as Sorbello 
did£ The choice Zi = 1 represents a proton in a jellium. The width of the square well potential was chosen to be 
equal to the screening length 1/ A and twice as large. The corresponding well depth was limited by the condition that 
just no bound state could be formed. The value of A is determined by the Fermi energy. While Sorbello chose five 
values for the Fermi energy, typical for metals ranging from sodium to aluminum, we have done the calculation for a 
whole range of Fermi energies. The results are plotted as a function of the Fermi wave number kp. The kp values of 
sodium and aluminum are indicated. 

Because A increases monotonically with the Fermi energy, the range of the corresponding screened Coulomb potential 
decreases with increasing kp, which reduction in strength is seen clearly in the bold solid curve. We compared v^k 
for the three potentials and found a clear decrease for the Coulomb potential with increasing kp, and a rather flat 
behaviour for the square well potentials, the one with 2/ A being markedly stronger than the one with the smaller 
width. The screening to second order in the impurity potential appears to be not negligible, but on the average 
15 ± 10%. For comparison we mention, that the result for a very small a value shows a screening which is at most 2% 
larger then the present result for a = r _1 = 0.01^ As a guide for the eye we gave the average of Z scr for the three 
potentials as a dotted line. Interestingly, this result does not imply that Sham's conclusion of a negligible screening 
is entirely wrong. If fact, he compared Z scr with Z W i n d- The ratio Z scr / Z W i n d is small indeed, but this comes from 
the large value of r. It may be clear that a comparison with Zi would have been more appropriate. 

In conclusion, we have shown that the starting formula of Bosvieux and Friedel gives the right driving force on an 
ion in a metal under the influence of an applied voltage. Because it always has been recognized as the first time that a 
quantummcchanical form for that force was written down, their contribution can still be characterized as a pioneering 
one. On the other hand, their prediction of a complete cancellation of the direct force has been falsified, because it 
was based on surface integral terms, which are zero due to the hermiticity of the system Hamiltonian. Further, the 
expression for the magnitude of the screening due to Sham does not give a negligible screening, but a screening of 
about 15%. This is in agreement with an earlier result based on another approach^ 

Taking all this together, the controversy regarding the magnitude of the direct force can be regarded as being 
resolved by now, and a unification of the various descriptions has been achieved. 
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